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Abstract 



Using perturbation theory, we explore the universal high momentum behavior of cor- 
relation functions of gauge invariant operators in planar noncommutative gauge theories. 
We find that the correlation functions are strongly enhanced when pairs of momenta 
become antiparallel. In particular, there is a transition from the previously noted ex- 
ponential suppression of correlation functions at high momenta to a more field theoretic 
behavior when the momenta of pairs of operators antialign within a critical angle. Some 
of our calculations can be extrapolated to strong coupling, and in particular we are able 
to reproduce precisely the supergravity prediction for the behavior of two point functions, 
including the coupling dependence. 



1 Introduction 



Many of the recent interesting developments in string theory have centered around various 
non-gravitational theories that arise from string theory in "decoupling" limits. These 
include local quantum field theories, such as the Af = 4 SYM theory, plus non-local 
theories without a known explicit description, such as little string theories. Somewhere in 
between are the non-commutative field theories which arise in certain decoupling limits 
involving D-branes with NS-NS two-form fields turned on in their worldvolume directions 
[1, 2, 3, 4]. These theories are non-local, yet possess explicit Lagrangian descriptions in 
the language of ordinary quantum field theory that allow perturbative calculations. 

Perturbative analysis of noncommutative field theories has revealed many interesting 
properties (see e.g. [5, 6]). In particular, the theories display a mixing of the infrared 
and the ultraviolet: short distance physics can lead to significant long range effects. 
From a technical point of view, these effects arise from internal momentum dependent 
phase factors which appear in non-planar graphs of the theory. On the other hand, 
planar graphs behave identically to their counterparts in the commutative theory apart 
from overall external momentum dependent phase factors. One might then expect that 
large N noncommutative field theories should have basically the same behavior as their 
commutative counterparts, since non-planar diagrams are suppressed by powers of 
However, it turns out that at least for noncommutative gauge theories, the commutative 
and noncommutative versions have very different behavior at large momenta, even in the 
planar, large N limit. 

This differing behavior can be seen easily for the particular case of M = 4 SYM 
theory by comparing the conjectured supergravity duals for the two cases. In the commu- 
tative case, the corresponding gravitational theory is of course type IIB string theory on 
AdS 5 x S* 5 which arises as the near horizon geometry of a stack of D3-branes [7]. In the 
noncommutative case, the dual theory is type IIB string theory in a background which 
has been worked out in [8, 9]. The two solutions have identical behavior in the "infrared" 
part of the space (near r = 0), suggesting that the two field theories are equivalent at low 
energies. On the other hand, the supergravity dual of the noncommutative theory has 
an Einstein frame metric [10] which is asymptotically flat in the r — > oo region (though 
the coupling goes to zero here) in sharp contrast with the asymptotically AdS dual of the 
commutative theory. Thus, we expect that the UV behavior of noncommutative M = 4 
SYM theory differs significantly from the behavior of its commutative counterpart, even 
for large N. The different behavior in the ultraviolet reflects the fact that unlike conven- 
tional field theories, non-commuatative gauge theories are not defined by a UV fixed point 
[11]. It seems particularly interesting to understand this UV behavior, since it provides 
the field theory dual of gravity on a space which is not asymptotically AdS. 

In this paper, we study the high momentum behavior of correlation functions of gauge 
invariant operators in noncommutative gauge theory. As discussed in [12]-[18] and re- 
viewed in section 2 of this paper, these operators involve open Wilson lines which extend in 
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the directions perpendicular to the operator momentum with a length proportional to the 
momentum. This transverse spreading at high momenta is in accord with the uncertainty 
relation associated with the commutation relations for the coordinates, [x M ,a^] = id^ v . 
The extended nature of the gauge invariant operators leads to qualitatively different be- 
havior of high momentum correlation functions as compared with the commutative theory. 

In [16], the two point functions of these operators were calculated at weak coupling 
and shown to have a universal exponential behavior with high momentum, 



where the dots indicate a regularization dependent term that cancels out in properly nor- 
malized correlation functions. The authors showed that the momentum dependence of 
this expression is identical to that of the two point function calculated using the super- 
gravity dual of noncommutative Af = 4 SYM theory, though the supergravity result has a 
different coupling dependence, obtained by the replacement g 2 N — > \ / g' 2 N in (1). In this 
paper, we show that even this coupling dependence can be reproduced by a perturbative 
field theory calculation extrapolated to strong coupling. We describe this result in section 
4 and offer two possible explanations for the agreement. 

In the rest of the paper we consider the high momentum behavior of more general 
correlation functions. We find a rather universal effect: the correlation functions are 
very sharply enhanced when the momenta of pairs of operators becomes antiparallel. 1 In 
particular, for a pair of operators with large nearly antiparallel large momenta of order k, 
(\k\ \6k\ ^> 1) we find that the correlation functions vanish rapidly when the angle between 
the momenta rises above some critical value, > p^j- As a result of this enhancement 
amplitudes fall off exponentially with the momenta for generic angles, but cross-over to 
a more field theoretic behavior below this critical small angle. 

The structure of the paper is as follows. In section 2, we review the construction of 
gauge invariant operators in noncommutative gauge theories and discuss in particular the 
operators corresponding to supergravity modes in the gauge theory - gravity correspon- 
dence for M = 4 SYM theory. In section 3, we provide a general formula for computing 
planar correlation functions of operators involving open Wilson lines. In section 4, we 
discuss the calculation of two-point functions and in section 5 we extend this discussion to 
higher point functions. We offer some comments and a summary of the results in section 



This is consistent with the interpretation of elementary quanta in the theory as dipoles, advanced in 
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[19]. 
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2 Gauge Invariant Operators 

In this section, we review the construction of gauge invariant operators in noncommutative 
gauge theories as well as the specific form of operators corresponding to bulk supergravity 
modes in the gauge theory - gravity correspondence for noncommutative j\f = 4 SYM 
theory. 

To begin, we recall the transformation properties of the gauge field in noncommutative 
gauge theory, 

-> U{x) -k -k U( X y - l -U{x) * d^{x) (2) 

We will take a U (N) gauge group, so that A^ is an N x N hermitian matrix. The gauge 
covariant field strength is then given by 

F^ u = d^A u -d u A ll + ig{A^A v -A u *A ll ) (3) 

and the (Euclidean) gauge invariant action for the gauge fields is 

S=^Jd 4 xTr(F flu *F flu ) (4) 

( equivalent ly, the star product may be replaced by an ordinary product). In more general 
theories such as the Af = 4 SYM theory, we may have additional matter fields, but the 
generalization of the action to a noncommutative theory is always obtained by replacing 
products by star products. 

We now turn to the construction of gauge invariant operators. A remarkable property 
of noncommutative gauge theories (with all fields transforming in the adjoint) is that 
translations in the noncommutative directions are gauge transformations, 

f(x" + e^K) = e ik ' x * fix 11 ) * e~ ik - x (5) 

As a result, it is impossible to construct local gauge invariant operators. On the other 
hand, it is possible to construct gauge invariant operators that are local in momentum 
space, as shown in [12, 16, 17, 18]. 

The construction employs open Wilson lines, defined by 

W(x,C) = Tr (p* e i9 fc A » dx ^ 

oo „! r i 

= YtiwT / da,--- da n (^---(^ 

Tr (A^(x + C(<7i)) *---*A IMn (x + CK))) 

(6) 

where P* denotes path ordering and C is a path parameterized by ((a) with £(0) = x. 
These transform covariantly under gauge transformations, 

W(x, C) -> U{x) * W(x, C) *U\x + C(l)) (7) 
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Figure 1: A general Wilson line operator with momentum k. 

Now, given a set of gauge covariant local operators {Oi} transforming in the adjoint, we 
may construct a gauge invariant operator local in momentum space by taking 

6(k) = J e ikx Tr (O^x) * W(x, d) * 2 (x + h) * • • • * W(x + h + . . . + Z n _i, C n )) (8) 

where U = 0(1)- Using the relation (5), it is simple to show that this operator is gauge 
invariant provided that 

+ . . . + % = e^K (9) 

This operator is precisely the result of taking insertions of any set of local covariant 
operators at arbitrary point on an open Wilson line whose endpoints are separated by 
k = 9^ u k u , as diagrammed in figure 1. 

In fact, it is somewhat misleading to refer to these objects as open Wilson lines, as 
the endpoints of the interval are not distinguished. Noting that 

J e ihx Tr (0 1 (x)*W(x,C)*0 2 (x + k)) = J e ik x Tr (0 2 (x)*0 1 (x)*W(x,C)) , (10) 

we see that the operators at each end of the open Wilson line may equivalently be taken 
at a single point by a cyclic rearrangement. 

For the case of noncommutative M = 4 SYM theory, there is now compelling evidence 
[16, 22, 23] that the operators corresponding to particle states in the dual gravitational the- 
ory are gauge invariant operators whose form is a special case of the construction described 
so far. Recall that for the usual commutative Af = 4 theory, operators corresponding to 
particles in AdS 5 x S 5 are chiral operators which take the form of a symmetrized trace 
of a product of gauge covariant objects (scalars, fermions, field strengths, or covariant 
derivatives of these), 

0(k) = J e ikx STr (B l ---B n )(x) (11) 

It turns out that the appropriate generalization of this operator to the noncommutative 
theory is given by 

0{k) = J e tk - x STy(B 1 * • • • * B n * W(x, k)) (12) 
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where W{x, k) is a Wilson line with a straight line path from x to x + k and STr denotes 
that the expression is to be averaged over all ways of inserting the Bs into separate points 
on the the Wilson line. Thus, the individual factors in the product become spread out 
linearly over a transverse distance k. 

In the following sections we study general correlation functions of these straight Wilson 
line operators, in particular looking for universal behavior at high momenta that arises 
from the extended nature of the operators. It is interesting to note that due to the 
asymptotically flat nature of the supergravity dual (in the maximally supersymmetric 
case), there are scattering states in the dual geometry, and it is likely that the correlations 
functions we compute correspond to S-matrix elements of these scattering states. This 
is similar to the situation in linear dilaton backgrounds [24], and is in contrast to the 
asymptotically AdS geometries. 



3 Correlation Functions of Gauge Invariant Opera- 
tors 

We have seen that the gauge invariant operators in noncommutative gauge theories have 
an extended structure that differs significantly from the corresponding operators in the 
commutative theory at high momenta. In this section, we investigate the effects of this 
extended structure on general correlation functions. 

Since we are mainly interested in the effects of the Wilson line part of the operators, 
we focus on the simple case of Wilson lines with a single insertion of a gauge covariant 
operator O. We define 

W(k) = J d 4 x e ik ' x * Tr (0{x) * P + e iff / A » dx ") 

= jr(tg) n ( d 4 x e ik ' x ( C M1 K) • ' -C^K) 

J J0<ai<-<a n <l 

Tr (0(x) * A^(x + CM) * • • • * A, n (x + CK))) (13) 

As before, P* denotes noncommutative path ordering, and the Wilson line runs over a 
path x^ = C M (°") , w hh C M (1) — C M (0) — ^ = <d^ u k u . We will mainly consider the case of 
a straight Wilson line for which C M (o") = k^a. 

We focus on the large iV theory, for which planar diagrams are dominant. For such 
diagrams, the star products appearing in the expansion (13) only have the effect of pro- 
viding an overall phase factor for the correlation function, which depends on the ordering 
of the external momenta. 

To focus on the UV behavior of the theory, we assume that the external momenta are 
asymptotically large. More precisely, we assume that the dimensionless quantity \k\\k\ 
is large. The amplitudes are then dominated at each order in the 't Hooft coupling by 
ladder diagrams, as discussed in [16]. These are diagrams, such as the one shown in figure 
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Figure 2: A typical planar ladder diagram in the correlation function of four Wilson line 
operators 

2, in which the gauge fields are contracted between various Wilson lines, with no internal 
vertices. 2 Since we also work at leading order in the 't Hooft expansion, we are further 
restricted to planar ladder diagrams. 

For planar diagrams in which the Wilson line gauge fields contract only with each 
other, it is straightforward write down a general expression for the correlation function. 
We find 

(w(h)---w(k n )) mia = s^^/^^w-m..)) 

Here, the first sum is over all ways of contracting the As between the various Wilson lines 
to give a planar diagram. The phase factor comes from the star products, but depends 
only on the ordering of external momenta since the diagrams are planar. In the second 
line, the integral runs over positions of the As on the various lines (i.e. the endpoints of 
the propagators), and the product is over all propagators in the diagram. In the integrand, 
x and x represent the two endpoints of a given propagator, while ( and ( are tangent 
vectors evaluated at x and x respectively. A typical diagram in the sum is shown in figure 
2. 

4 Two Point Functions 

We start by calculating two point functions of the gauge invariant operators (13). The 
calculation was performed previously for weak coupling in [16]. Here, we note that within 
the approximation of summing only over ladder diagrams, the result may also be eval- 
uated at strong coupling. We find that this strong coupling result exactly reproduces 

2 Briefly, the dominance of ladder diagrams at large momenta arises because the endpoints of each 
"rung" are integrated over the (long) Wilson lines, giving rise to factors of k. These factors enhance the 
ladder diagram relative to any other diagram of the same order in the 't Hooft coupling. 
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the behavior of the supergravity result including the coupling dependence. At the end of 
this section we offer two possible explanations for the agreement. A similar result was 
obtained in [20] in the context of Wilson line correlators for the commutative M = 4 
theory, see also [21]. 

For the two point functions, factoring out the 5-function of momentum conservation, 
the general formula (14) becomes 

W{k)W\k)) = f rf 4 Ae* A /(A)E(APr ft / dn I do; t T7~~~T7~ ^ ^ 

J n i=l |A T fe(T,: — <Ti)\ 



Here A = j^-, A = x\ — x 2 is the relative position of the operators, and we define 

(0 1 (x 1 )0 2 (x 2 )) = f(x 1 -x 2 ) . (16) 

Because of the restriction to planar ladder diagrams, the integrals over a and r (which 
describe the endpoints of the various rungs) are ordered, < a± < • • • < a n < 1, 
< T\ < ■ ■ ■ < r n < 1. In order to extract the universal large momentum behavior arising 
from the Wilson lines, we follow [16] and set /(A) = 1. 

The integrand \ A+ ^._ a .^2 is a steep function (for large |A;|) which is maximized at: 

= -pr (17) 

| A; I 

where Ay is the component of A parallel to k. In order to attain this maximum within 
the integration region one requires: 

< A), < \k\ (18) 

Since the integrand is steep, we would like to evaluate it using a saddle point ap- 
proximation starting with the Tj integrals. However, due to the ordering of the points, 
T i < T i+i, the peak of the integrand may be close to the edges of the integration regime, 
so the saddle point approximation is not automatically justified. The required condition 
is that the width of the integrand as a function of a given r iy 

IA, I 



1*1 

is less than the average spacing between the r^'s, K Thus, to use the saddle point 
approximation, we need 

i^^Ul (19) 



When this condition is satisfied, defining 

AA; 

fJ[J "J '\A + k(Ti-a. 



4=n/d Ti /^ (20) 
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we find (first evaluating the r integrals using saddle point and then performing the a 
integrals directly), 

l n „ \ * (21) 

To see that this estimate cannot hold when 77 > 1 , note that the integral I n is bounded 
from above by the maximum value of the integrand times the volume of the integration 
region 

In ~ (dp A^ (22) 

It is easy to see that the saddle point result (21) exceeds this bound when rj > 1. In fact, 
the upper bound (22) provides a good approximation to I n for rj > 1 (we will show this 
explicitly in section 5 using a lower bound on /„). 

In either regime, the integral over Ay (which we restrict to the range (18)) can be 
readily performed, remembering that k\\ = 0. Suppressing the subscript of Aj_, we get: 

W{k)W\k)) ~ \k\ J d 3 Ae* A ]^(^-) n / n (23) 

{>lL (k\ n nA <r 1 

"! VA^ k ^ f OA \ 

_A^_ (R\ 2n nA . 1 ( 24 ) 

(n!) 2 VAj fc > 1 

The correct estimate for I n depends now on the values of A, n dominating the ex- 
pression for the two-point function. For both estimates of /„, the integration over A is 
dominated by the regime |A| ~ therefore one has 11 ~ in the relevant region of 

integration. Thus, for n < \J~kk we may use the 77 < 1 estimate when integrating over A, 

while for n > \J~kk, we should use the 77 > 1 estimate. 

As explained in [16], we should extract the contribution which is non-analytic in k 
for each n, since the analytic parts correspond to contact terms in position space. The 
results for the relevant Fourier transforms can be found in [16]. We are left with a sum 
of terms J2 n Jn where the n dependence of J n is given by 

\ n (kk) n 



where 



i' " " ~ ft < C V kk 

(n!) - , >- 25 

(n!) 2 (2n)! n > V KK 

For each value of the parameters, the terms in this series will increase to a maximum 
value at some n = M and then decrease, with the terms for n 3> M giving negligible 

contribution. At weak coupling, we find that M <^ \f~kh so only the small n form of J n is 
relevant. Performing the sum using the small n estimate, we recover precisely the result 
(1) of [16]. 
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On the other hand, at strong coupling, M ^> y kk so the series is dominated by terms 
of the large n form. Summing the series in this case, we find (ignoring prefactors) 



We note that up to numerical factors (which we have not been careful about), this 
result is precisely equal to the supergravity result [16] including the coupling dependence. 

We offer two possible explanations for this agreement. Firstly, in the case of chiral 
operators for the Af = 4 theory, it is reasonable to guess that two point functions should 
obey a non-renormalization theorem such that the leading perturbative result is valid at all 
values of the coupling. But in our case, the leading perturbative result is precisely this sum 
over ladder diagrams. The higher order ladder diagrams are not perturbative corrections 
but actually the leading contribution to the correlation function from the higher order 
terms in the expansion of the operators (13). Thus, assuming a nonrenormalization 
theorem, we would expect that the sum over ladder diagrams at strong coupling should 
reproduce the supergravity answer. This is exactly what we find, so our result may be 
interpreted as providing evidence for this nonrenormalization theorem. 

A second possibility for the agreement is that as noted earlier, at high momenta the 
ladder diagrams give the leading contribution order by order in the 't Hooft coupling, 
while diagrams with internal vertices are suppressed by inverse powers of momenta. It 
is therefore possible that in the large momentum limit, the restriction to a sum over 
ladder diagrams is sensible even at strong coupling, providing another possible explanation 
of the agreement we find. If this is the case, it would indicate that both chiral and 
nonchiral operators in the theory share this universal behavior at large momenta. From 
the supergravity point of view, this would indicate a similar behavior in the UV region 
for both supergravity modes and stringy modes. 



We turn now to higher point functions of the noncommutative operators. We will show 
that for large momenta, the correlation functions are sharply peaked, preferring the kine- 
matics in which pairs of momenta are antiparallel. The properly normalized correlation 
functions generically exhibit exponential decay at large momenta, as described in [16]. 
However, at small enough angles we find that the amplitudes are enhanced, and the 
dependence on momenta has a cross-over to a more typical field theoretic behavior, as 
described below. 

We consider then the correlation functions (Wi(ki)W2(k 2 ) ■ ■ ■ W n (k n )). In the limit of 
large hi in the sense described above, we will show that ladder diagrams must generally 
be included even at weak coupling, as was the case for the two point functions. As we 
have discussed, the ladder diagrams dominate any diagrams with internal vertices at each 




(26) 



5 Higher Point Functions 
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order in the 't Hooft coupling. Diagrams with internal vertices are suppressed by powers 

of 777777, so working in the order of limits in which the momenta are taken large first, 
|fc||fc| 

we may consider only the ladder diagrams. These may be summed explicitly both for 
weak and strong coupling, and we provide the results for both cases. We emphasise that 
summing over ladder diagrams at strong coupling can be justified only with our carefully 
chosen order of limits. The agreement with the gravity results for the two-point function 
is an indication that the sum over ladder diagrams may be sensible even in the limit of 
strong coupling. 

In subsection 5.1 we derive the conditions under which it is necessary to include 
ladder diagrams connecting a given pair of operators W\, W 2 . We will see that in a weak 
coupling calculation and for generic angles between the momenta hi, there is no need to 
include ladder diagrams, and we may use the leading perturbative result. However, when 
two momenta ki,k 2 are nearly anti-parallel, the sum over ladders connecting Wi,W 2 is 
necessary, and provides an enhancement of the amplitude. At strong coupling, ladder 
diagrams will be important for all angles. 

In section 5.2 we present the approximation schemes which are appropriate for eval- 
uation of this universal behavior in both the weak and strong coupling regimes. These 
approximation schemes generalize the two estimates (77 < 1 and r\ > 1) for the two point 
functions, presented in the previous section. In the sections 5.3 and 5.4, we perform the 
calculation of the ladder diagram contributions in each of these regimes. 

The reader interested only in the results of the calculation is invited to skip to section 
6, where our results are summarized and discussed. 



5.1 Conditions for Including Ladder Diagrams 

Given a general correlation function, we would first like to determine the conditions under 
which it is necessary to include ladder diagrams connecting a pair of operators Wi(k) and 
W 2 {— I) in a perturbative calculation. 3 We will find that ladder diagrams dominate the 
leading order perturbative result when the angle between k and / is less than a particular 
upper bound. 

To begin, consider a diagram involving an additional contraction between gauge fields 
in the Wilson lines associated with W\ and W 2 . Relative to the leading order perturbative 
result (with no contractions between the Wilson lines) this diagram will have an additional 
factor of 

X(k ■ I) f da f dr 1 (27) 

Jo Jo (A + ka - It) 2 

We wish to determine when this factor is of order 1 or larger, in which case the diagram 
(and higher order ladder diagrams) must be included. 

3 We have written the momentum of W2 as — I because this will simplify the equations. 
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Figure 3: Angles and momenta for which ladder diagrams dominate the leading pertur- 
bative result. 

For simplicity, we assume that A is perpendicular to k and I, which are taken to 
have the same magnitude. In this case, the integral may be evaluated explicitly, and the 
relative factor becomes 

Acot^sinh' 1 (^^j , (28) 

where <f> denotes the angle between k and I, assumed to be small. The distance A lies 
purely in the commutative directions. 

The typical value of A for operators of momentum k will be roughly A = so we 
estimate that the ladder diagram contributions will be important for 



4sinh- 1 (|A;||fcj0) > 1 





(29) 



This region is diagrammed in figure 3. We see that ladder diagrams must be included 
for small angles whenever \k\\k\ > 1/A. In particular, for momenta such that \k\\k\ ^> 1/A, 
ladder diagrams must be included when 

^<Alog(|A#|). (30) 

At larger angles, the ladder diagrams do not make a significant contribution, and we can 
trust the leading order perturbative result. At weak coupling, this cutoff angle remains 
small (^C 1) as long as \k\\k\ < e*, so the leading order perturbative result is fine for 
generic angles (0 ~ 1). On the other hand, it is obvious that ladder diagrams will 
dominate the leading perturbative result at strong coupling for any angles. 



5.2 Approximation Schemes 

We would like to calculate the contribution to the correlation function 

{W 1 {k)W 2 {-l)---W n ) 



(31) 
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Figure 4: Simple configuration of Wilson lines. 

from the sum over ladders connecting the Wilson lines in W\, W 2 in the case where k and 
—I are very large and nearly antiparallel. We are interested in universal behavior arising 
from the extended nature of the Wilson lines, independent of the details of the operators 
W t . 

From the general formula (14), we see that the correlation function has the following 
dependence on k, I 

(W 1 (k)W 2 (-l) ■ ■ ■ W n ) ~ / d 4 A exp(*A • (k - Z)/2)/(A, • • •) E J ™ (32) 

where / represents the commutative position space correlator, and X n is the contribution 
coming from the ladder diagram with n rungs connecting W\ and W 2 (see figure 4), 

rl r<?2 rl rT2 _ _ *L A 

I n = da n --- da, / dr n ---\ dn (k ■ l) n JJ — = ^- (33) 

Jo Jo Jo Jo iJ- 1 \/^ + ka i -lT i \ 2 

We would like to determine the behavior of the integral Z n for large k and — / equal in 
magnitude (for simplicity) and nearly antiparallel. In the full expression for the momen- 
tum space correlator, the separation A is integrated over, so we are mainly interested in 
the behavior of the integral in the region of A where the integral is largest. 

Let A^ and Ay be the components of A perpendicular to and parallel to the plane of 
noncommutativity. 4 The largest values of X n occur if Ay is such that the two lines cross 
in a projection onto the plane of noncommutativity. From figure 5 it is easy to see that 
this region of Ay has area \k\\l\ sin(0), where <fi is the angle between k and /. 

For Ay in the region described, the minimum distance between the two lines will 
be |Aj_|. We expect that the main contribution to the integral will be from the region 
in which most of the propagators have length of the same order of magnitude as this 
minimum length. In this case, most of the propagator endpoints remain within a distance 

4 For unitarity, we assume that the noncommutativity parameter O'"' is nonzero only for a pair of 
spatial directions and refer to these directions as the plane of noncommutativity. 
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Figure 5: Maximally contributing configurations of two Wilson lines. 



of roughly '—^ from the crossing point. The average separation between endpoints on 
each line will therefore be roughly min ( ^jp, — ). 



If the average separation between the endpoints is much greater than the average 
propagator length, it is a good approximation to extend the region of integration for each 
r between — oo and oo. This approximation will therefore be valid if 



We will call this the small n approximation, since there is an upper bound on n (though 
it may be large). This generalizes the 77 < 1 regime in our evaluation of two point function 
for which the saddle point approximation was justified (21). For n larger than this bound, 
we will require other methods to evaluate the integral. 

We now turn to a detailed evaluation of the integral X n in both the small and large 
n approximations. It will turn out that for weak coupling calculations, only the small n 
regime is relevant, while for strong coupling, we require the large n behavior. 

5.3 Small n Approximation 

We begin by evaluating the integral X n in the small n regime. As discussed, for n satisfying 
the bound (34) it is a good approximation to extend the integration region for the r's to 
the whole real line for a typical set of cr's (equivalent to a saddle point approximation). 





(34) 
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This also provides an upper bound on the integral for all values of n and the other 
parameters. 

After performing this saddle point evaluation, the remaining integral takes the form 

1 n = \ da n ---\ dv 1 f(a 1 )---f(a n ) = -(/ daf(a)) n 
Jo Jo nl Jo 

= ^-H n (35) 
nl 



where 



m = / 



dr- 



\k ■ I 
|A + ka — lr\ 
Xk-l 



l 2 {A + ka) 2 -{l-{A + ka)f 



(36) 



We can perform the the integral over a exactly, though the result is somewhat messy. In 
order to get an idea of the dependence on and k, it is illuminating to consider the case 
in which the lines have equal length. 

The maximum contribution comes from the value of Ay for which the lines cross 
symmetrically in the projection to the noncommutative directions. This occurs for A- k — 
—k 2 sin 2 0/2, and in this case, the result is 



k sin <f> \ 2 ill fc sin 4> 

iJ = Acot01og| Y^g^ ?J (37) 



2A ) ' 2A 

Here, once again for the sake of brevity we denote by A the distance in the commutative 
directions, previously denoted as Aj_. Though we have chosen a particular value for Ay, 
we note that X n should be roughly constant over the region of Am shown in figure 5 since 
the dominant configurations in the integral will have all propagators clustered around the 
crossing point. 

Recalling that T n = H n /n\, we see that for angles j, 

^ 4 (4' 1 -I£+->)" < 38 » 



while for 3> —, we get 



*-a( A H$) (39) 



These formulae accurately describe the behavior of X n in the small n regime, (34). 
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5.4 Large n Approximation 



We now consider the case in which the small n approximation is no longer valid. This 
occurs when the average separation between points on each line becomes less than the 
typical distance between the lines. In this limit, the propagator endpoints become dense 
on the two lines, and we expect that most of the contribution to the integral comes from 
a small region about some preferred configuration. 

To get an idea of the k and <fi dependence of X n in the large n regime, we again take the 
simple configuration with \k\ = \l\ and Ay = (this is the configuration shown in figure 4 
if the separation A in the diagram is understood to be perpendicular to the lines). Since 
the integrals are difficult to evaluate directly in the large n regime, we will determine the 
behavior of I n by establishing upper and lower bounds. 



5.4.1 Large n, Upper Bound 

First, we will determine an upper bound on X n for this configuration using the relation 

j\fl-f2---fn\<(j\fl\ n )---(j\fn\ n )" (40) 

where /$ are arbitrary functions on some integration region (this is a generalization of the 
Cauchy-Schwarz inequality). This implies that 

Tn<(f[T m )" (41) 
m=l 

where X m is obtained from X n by replacing the product of propagators in the integrand 
by n copies of the mth propagator. The integrals T m are much easier to evaluate, as we 
will see. 

The integrand of X rn depends only on a rn and r m , and it is straightforward to perform 
the integrals over all other as and rs. We find 

± m = (X~k ■ l) n C da m f 1 dr m , A /( y m)/( ^ m) , 2 (42) 
A) Jo \A + ka m -lr m \ 2n 

where 

j.m— i\n-m 

'<«> = (m -l)Kn- m )l (43) 

We note that f(t) is sharply peaked at t — while the denominator of T m is sharply 
peaked near a = t = 0. The integral will receive most of its contribution from these two 
regions. For small enough A, it is the latter region near the origin that dominates, since, 
as we will see, X n diverges as A — > 0. It is convenient to rewrite the integral in terms of 
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cr + = a + t and er_ = o" — r. This gives 

~ A n fc 2ri cos" 

m ~ [(m-l)!(n-m)!] 2 { ' 

1 ^ ^ (1 ((T 2 _ ^))m-l ((1 _ 1^)2 _ l a 2_ )n -m 



2 Jo " (A 2 + Fcos 2 (fy 2 +Psin 2 (f)a 2 )™ 

(45) 

The <t_ integral is sharply peaked at o"_ = and may be evaluated using the saddle point 
approximation. Near a + = 0, the result (for m > 1) is 



f A n P"cos n /" j_ <r+ m_1 (l - ^^ 2 "" 2m 



-lJo + (A 2 + k 2 sm 2 (£)al) n 



2 2m-l[( m _ l)!( n _ m )!]2 y m _ I J Q f ( A 2 + fc 2 sin 2 (§)(7 2 )» 



The dominant term for small A (A < k(f>) may now be computed exactly, and we find 
(for 1 < m < n) 

m [(n-l)\] 2 \m-ljA m (n-m)V^ zr Tsm 2m ^ \A 2 J 1 } 

Finally, using the relation (41) we find an upper bound for > k/A (ignoring constants) 



*■ £ W l^fj (48) 

Proceeding from equation (46) for < kj A we find 
5.4.2 Large n, Lower Bound 

We now determine a lower bound on X n in the large n regime. To do this, we restrict the 
integration region by choosing an ordered set = a < a 1 < • • • < a n — 1 < a n = 1 such 
that 

at < ai, n < a i+ i (50) 

Then certainly, 

fo-i+i ra.i+1 \k ■ I 
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This bound is valid for all choices of a i: so in order to establish the best lower bound, 
we would like to maximize the last expression over all such choices. 5 Despite a great 
reduction in the size of the integration region, this approach actually gives a reasonable 
estimate of X n for large n since the full integral is dominated by contributions from a very 
small region about some preferred configuration, and the maximization over choices of a« 
essentially selects this dominant region. 

Since we take n to be very large, the values a; are very closely spaced. We therefore 
can define a continuous functions f(x), such that: 

= /(-) (52) 

Th 

This satisfies the boundary conditions: 

/(0)=0 /(l) = 1 (53) 

Taking logarithms of both sides of (51) and replacing the sum by an integral we then 
obtain 

,og a gjj) • (54> 

where 

G n {f) = n f 1 dx{2\og(f'(x)-)-\og(\A + kf(x)-lf(x)\ 2 )} (55) 

J Th 

= -2n\og(n)+n ^ dx{2 hg(f'(x) - log(A ± 2 + f 2 S 2 )} (56) 
Jo 

Here, we have defined S = 2\k\ sin(|), and to avoid clutter, we omit the subscript in Aj_. 
We now maximize G by solving the Euler-Lagrange equation for /. This gives 

21og(/')-log(A 2 + / 2 S 2 ) = C (57) 

where C is a constant. 

Solving this and choosing integration constants to satisfy the boundary conditions 
(53), we find 6 

A S 
fmin(x) = — sinh(a: • sinh _1 (— )) (58) 

It is now simple to plug this value of f min into G to obtain a lower bound 

I n > (Xk ■ f)" e G(/mi " ) 

A n cos n / 1 . , 1 (2ksm^\\ 2n 



5 Note that we have restricted the a and r integrals in the same way because we are still assuming a 
symmetric configuration in which k and I have the same length with A|| = 0. In more general cases, it 
would be necessary to divide the a and t integrals differently to obtain the best bound. 

6 Note: sinh _1 (a;) = log(x + \f\ + x 2 ) 
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For k<p <C A , which includes the exactly antiparallel case, this gives 

X„^(|)"( 1 -i(^ + ...| (60, 

while for k<p 3> A (the non-parallel case for small enough A) we get 

A n cos"0 / 1 , (4~ksm£\\ 2n 

^^(s^H^JJ (61) 

Note also that 

n 2n = e 2(nln(n)-n) e 2n _ (q 2 ) 

so for example in the antiparallel case, we have upper and lower bounds 

11 11 

<^<7-tt^ (63) 



(n!) 2 (eA) 2n ~ " ~ (n!) 2 A 2n 

This justifies the claim in section 2 that the upper bound provides a good estimate of 
I n in the large n regime for the antiparallel case. 

We have now determined the behavior of X n for general values of n and the other 
parameters. In the next section we summarize these results and discuss the physical con- 
sequences for the behavior of the correlation functions at various values of the parameters. 



6 Summary of Results 

We have considered correlation functions of the form 

(W^k^i-V-'-Wnikn)) (64) 

and noted that even at weak coupling, ladder diagrams connecting W 1 and W 2 must be 
included when the momenta k and —I are nearly antiparallel, as summarized in figure 3. 
For the ladder diagram contributions, the dependence on k and / is given by 

(W 1 (k)W 2 (-l) ■ ■ ■ W n ) ~ / d 4 A exp(*A • (k - Z)/2)/(A, • • •) E J ™ (65) 

J n 

where / represents the commutative position space correlator and X n is defined in (33). 

Our results for the behavior of X n are summarized in figure 6. For momenta of magni- 
tude k with a given fixed angle between the, we find four regions of behavior as n and 
A are varied: 



18 



n 





larg 


e n regime 






c 




D 




/a 

smal 


n regime 


B 




1 



Figure 6: Regimes of behavior for X n in the A — n plane. 

small n: 



In this regime, defined by 

1 k 

n < min(-, — ) , (66) 

the typical spacing between points on the Wilson lines is greater than the average sepa- 
ration between the lines (propagator length). We find that 

T - J n - V A V 24A2 + ^ A (j> f R7 \ 

" U (****<¥))" 4>* ( ) 

These correspond to regions A and -B respectively in figure 6. 
large n: 

In this regime, for which 

1 k 

ra > min(-, — ) , (68) 

the points are densely packed on the two Wilson lines. For ^ > ^ (region D) we have 
determined upper and lower bounds given by 

A n cos n / 1 , /4fcsin|\\ 2n ^ A" f kcos<p\ n , . 

^(^ io H^JJ - z " s wu4) (69) 

Finally, in region C with ^ < ^ we find 
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In evaluating the momentum space correlators, we must sum over n and integrate 
over A as in equation (32). In this case, because of the additional A dependence of the 
function /(A) in (32), it is simplest to consider first the sum over n with everything else 
fixed. 

Note that as a function of n, the series of terms I n is exponential (X n ~ in the 
small n regime and then falls of faster (J n ~ 7^2 ) m the large n regime. Thus, it is obvious 
that the series converges for any values of the parameters. Furthermore, just as for the 
two point function, there will be some M such that the terms start decreasing for n > M, 
and terms with n ^> M will give negligible contribution. 

At weak coupling (A 1), we find that M lies well within the small n regime as long 



7 

as 



<P<te^ . (70) 
k 

In this case, the large n terms have negligible contribution and we find 

£2,2 

£J n = e A A. e 53A7 < J 

£Zn=(?)* f «0<Alog| 

For larger angles, ladder diagrams are no longer important and the result will be given 
by the leading order in perturbation theory. 

For strong coupling A 1, M lies in the large n regime. In this case, the series is 
dominated by terms in the large n region, and we find (for small angles) 

e# e _l (-) <EI n <efe 4 W 0« A 



(*) f <EX n <exp(^) 0»f 

The detailed form of the correlation functions would now be obtained by inserting 
these expressions for T n into the expression (32) and integrating over A. However, the 
universal behavior at large momenta, including the angular dependence, is already quite 
clear from the various expressions for X n . 

We see that both for the weak coupling and the strong coupling expressions, £X n is 
very sharply peaked near = 0, and decays very rapidly for > A. Roughly, the typical 
value of A in the Fourier transform will be 7 so we see that there is a critical angle = tt 

k & ^ kk 

below which the correlation functions are strongly enhanced. 8 For = we get precisely 
the exponential dependence on momenta that was observed in the two point functions. 

7 For ^ < ei, this will be true at all angles at which ladder diagrams are important. 
8 We should point out that the external momentum dependent phase factor in (14) gives rise to an 
additional angular dependence, cos(kksin(J)) which oscillates rapidly for large momenta. 
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In [16], it was pointed out that in the absence of a natural normalization of individual 
operators the most meaningful objects to consider are normalization independent ratios 
of correlation functions, 

(W 1 (k 1 )W 2 (k 2 )---W n (k n )) 
(W 1 (k 1 )W 1 (-k 1 )) 1 2...(W n (k n )W n (-k n )} 1 2 

They argued that the exponential behavior of the two-point functions would generally 
lead to exponential suppression of the correlation functions at high momenta. Given our 
results, we see that this will be true for generic angles, however, if the momenta of two 
operators become pairwise antiparallel (0 i), the numerator will also exhibit exponen- 
tial behavior that should precisely cancel the exponential behavior in the denominator. 
In this case, the momentum dependence is no longer universally determined by the Wil- 
son line contribution. Rather it is determined by the correlation functions of the local 
operators attached to the Wilson line. Therefore we expect that it reverts to more typical 
field theoretic behavior. In particular, in the planar limit the behavior is determined by 
the UV fixed point of the commutative version of the theory. It would be very interesting 
to understand the meaning of the angular dependence we have observed from the point 
of view of the gravity dual, but we leave this as a question for future work. 
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